Effects of Scale-Free Disorder on the Anderson Metal-Insulator Transition 
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We investigate the three-dimensional Anderson model of localization via a modified transfer- 
matrix method in the presence of scale-free diagonal disorder characterized by a disorder correlation 
function g{r) decaying asymptotically as r~". We study the dependence of the localization-length 
exponent f on the correlation-strength exponent a. For fixed disorder W, there is a critical Oc, such 
that for a < Oc, v = 2/a and for a > a^, v remains that of the uncorrelated system in accordance 
with the extended Harris criterion. At the band center, v is independent of q but equal to that of 
the uncorrelated system. The physical mechanisms leading to this different behavior are discussed. 
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The successful analysis of the metal-insulator transi- 
tion (MIT) in the Anderson model of localization jlj 
has hitherto been limited to short-range or uncorrelated 
diagonal disorder. In this Letter we report the effects 
of long-range power-law correlated disorder — so called 
scale-free disorder — on the MIT. Scale-free disorder is 
omnipresent in nature. It is found in many diverse situa- 
tions in biological 0,0 and in physical H|p systems, in 
city growth patterns [6| and in economics [TJ . The effects 
of scale-free disorder on the critical properties of physical 
systems have recently received much renewed attention 
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For the Anderson model of localization previous inves- 
tigations of scale-free disorder have concentrated on the 
one- (ID) [il [H 113 and two-dimensional (2D) jH El 
cases. For the ID Anderson model, it has been shown 
that for energies close to the band edge the presence of 
scale-free diagonal disorder causes states to be strongly 
localized ^3J, while at the band center the states tend 
to have localization lengths ^ comparable to the system 
size the 2D Anderson model with scale- free dis- 

order an MIT of the Kosterlitz-Thouless transition type 
[21I [2^ has been observed. 

Of particular interest to us is the influence of scale- 
free disorder in the neighborhood of an MIT when the 
localization length ^ becomes sufficiently large. Scale- 
free disorder can affect the character of the divergence 
as shown in Refs. 0, for the classical percolation 
problem. The principle finding of Refs. |24i |25j is that 
the critical exponent v, governing the divergence of ^, in 
a scale-free disordered classical system can change when 
the correlator g(r — r') falls off with distance as a power 
law, i.e., cx |r — r'|^". The critical exponent of the clas- 
sical percolation u = 4/3 crosses over to ly = 2/ a for 
a < 3/2, i.e., when the decay of the correlator is slow 
enough. 

The Harris criterion l2d l27| summaries the effects of 



short-range correlated disorder on a critical point. The 
criterion states that v for the disordered system and the 
clean system are identical provided that dv — 2 > 0, 
with d being the dimensionality of the system. The in- 
equality is derived by demanding that the fluctuations 
of the random potential within a volume given by f do 
not grow faster than their mean value as the transition 
is approached. For power-law correlated potentials, like 
the above scenario, there is an extension of this criterion 
originally suggested in Ref. and further studied in 
Refs. |l6l |25| . The behavior of v in the presence of scale- 
free disorder is well described by the extended Harris 
criterion which can be stated formally as 



2/q; if a < ac 
if q; > die 
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where is the critical exponent without correlations in 
the disorder. Eq. implies that there is a well-defined 
critical value ac = 2/1^0, below which correlations are rel- 
evant and above which correlations are irrelevant. Nu- 
merical studies for 2D classical percolation pj|] are indeed 
in good agreement with Eq. Quite recently, numeri- 
cal investigations of long-range correlations in models of 
2D quantum-Hall systems corroborated its validity for 
the quantum case 0, 0| . Therefore it appears to be 
possible that the criterion based on potential fluctuations 
can be applied also to quantum phase transitions such as 
the Anderson-type MIT, where v is determined from the 
divergence of the quantum-locdXvL&tioTi length at the crit- 
ical point. It is the purpose of this Letter to investigate 
this possibility. 

We define ve and uw as critical exponents of the local- 
ization length ^ oc |i? — i^d^"^ at fixed disorder strength 
W and ^ oc \W — Wc]"""' at fixed energy i?, respectively. 
We use the symbol v to denote both exponents ve and 
vw- For fixed W and a < ac, the critical exponent ve 
obeys the extended Harris criterion, as shown in Fig. ^ 
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binding Hamiltonian 28] in site representation 



FIG. 1: The localization-length exponent ue{q) as a function 
of the correlation-strength exponent a a,tW — 12. Error bars 
reflect one standard deviation. The horizontal line indicates 
the uncorrelated vo = 1.66 ± 0.06, the vertical dotted line is 
ac = 1.21. The dashed line for a < Uc gives the extended 
Harris criterion Q. The grey areas denote error bounds of 
one confidence interval arising from the error in vq. Deviations 
from the extended Harris criterion for small a < 0.5 are due 
to finite-size effects. 
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FIG. 2: The localization-length exponent uw{oi) as a function 
of the correlation-strength exponent q at iJ = 0. Error bars 
reflect one standard deviation. The horizontal line indicates 
the uncorrelated vo = 1.49 ± 0.03, the vertical dotted line is 
ac = 1.34. The dashed line for a < Qc gives the extended 
Harris criterion Q. The grey areas denote error bounds of 
one confidence interval arising from the error in j/q. 



At the band center E = Q, however, u-^^ remains inde- 
pendent of a as can be seen in Fig. |21 This means that 
scale-free disorder increases the critical exponent ve f^or 
small a while leaving unchanged. 

Our calculation is based upon the Anderson tight- 
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where (i, j) denotes a sum over nearest-neighbors and \i) 
is an atomic- like orbital at site i. The random on-site po- 
tentials Si are chosen from a Gaussian distribution with 
zero mean and variance W"^ 112. For uncorrelated Gaus- 
sian disorder, the dependence of Wc[E) is known [2^ . 
in particular Wc(0) = 20.9 ± 0.5. We generate scale- free 
disorder by use of the modified Fourier-filtering method 
(FFM) as outlined in Refs. 113,1^, so that the random 
on-site energies have an asymptotic correlation function 
(eiEi+r) ^ f " hi real space. The average is done over 
spatial positions and many disorder realizations. We note 
that large a corresponds to the nearly uncorrelated case 
which shall serve as our point of reference, while small a 
is the strongly correlated case. 

In principle the usual iterative transfer-matrix method 
(TMM) Hi m H S m allows us to determine the 
localization length A of electronic states in a quasi-lD 
system with cross section M x M and length L ^ M, 
where typically a few million sites are needed for L to 
achieve a reasonable accuracy for A. However, the use of 
the non-iterative FFM procedure to generate scale-free 
disorder necessitates a complete storage of the on-site 
potentials and consequently, the iterative advantage of 
TMM is lost as computer memory requirements become 
rather large. 

In order to circumvent this problem, we have modified 
the conventional TMM. We now perform the TMM on a 
system of fixed length Lq of the quasi-lD bar. After the 
usual forward calculation with a global transfer matrix 
, we add a backward calculation with transfer matrix 
72"^ . This forward-backward- multiplication procedure is 
repeated K times. The effective total number of TMM 
multiplications is L = 2KLq and the global transfer- 
matrix tl is 



TL - (^^.•T^••T^„Ti„.•.^„.•.Tl) 



K 



K 



As usual, we diagonalize the matrix 

l/AKLo 



lim 
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This modified TMM has previously been used for two 
interacting particles j37i|. 

After establishing Eq. the calculation of A follows 
that for the conventional TMM [38]. The matrix 
is symplectic with paired eigenvalues exp(±7) with 
Lyapunov exponents 7 . Physically, 7 determines the in- 
crease or decrease of the envelope of the wave function at 
long distances. The localization length is defined as the 
inverse of the smallest Lyapunov exponent, A = l/7min 
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[3^ liol liH . For a reliable convergence check, only the 
accumulated changes to A after each complete forward- 
and-backward loop need to be taken into account, not 
all such changes in the bulk (n = 2, 3, . . . , — 1) of the 
sample. Usually, convergence can be achieved after just 
a few K multiplications. 

The critical behavior can be determined by numerically 
establishing for the reduced-localization lengths A = 
A/M the one-parameter scaling hypothesis 
as A(a;,a) = T[£,{x,a)/M] in the vicinity of the MIT, 
where ^(a;, a) is the three-dimensional (3D) localization 
length in the thermodynamic limit. And accordingly 
[ll l34j . ^ diverges when the tuning parameter x, which 
can be the disorder strength W or the energy is tuned 
to its critical value Xc, i.e., Wc or Ec as 



^(x) ~ |a; - Xc[ 



(5) 



The value of v is estimated from the one-parameter 
scal ing hypothesis 0, 01 by finite-size scahng (FSS) 
[i^ liil l44l |. The FSS procedure performed here follows 
closely the approach in Refs. jil, jii, [ij [i^. The basic 
idea [33 is to construct a family of fit functions which 
include corrections to scaling [Sg due to an irrelevant 
scaling variable and due to non-linearities of the disorder 
dependence of the scaling variables ■ 

Numerical results have been obtained for samples of 
lengths up to at least Lq = 1000, system widths M — 
5,7,9,11,13 and for various a values. For band cen- 
ter data as in Fig. |21 widths up to M — 15 have been 
used. For each sample, convergence was assumed when 
the value of 7 changed by less than 10~^ after a complete 
forward-and-backward loop. We emphasize that since no 
self-averaging is used, the more stringent standard TMM- 
convergence criterion (3^ is unnecessary. With at least 
100 samples for each A(W, E, a), we find a relative error 
of the sample-averaged localization lengths of < 5%. 

The phase diagram of localization. Fig. O shows the 
mobility-edge trajectories, which separate extended from 
localized states, for different values of a. The symmetry 
with respect to ii^ = holds as in the uncorrelated case. 
Wc(a) at ii^ = increases monotonically as a — > 0. In the 
absence of correlations, states in region I are extended 
whereas states in II are localized. The phase bound- 
aries between the metallic states (region I) and the insu- 
lating states (region II) are modified in the presence of 
scale-free disorder. With uncorrelated disorder, states in 
region III are localized by quantum interference effects 
[23,|53|- Correlation in the disorder potential will lead to 
a smoother disorder and a decrease of interference and 
thus increased localization lengths. Indeed, we find that 
for a < cxD, states in region III become extended and the 
phase boundary is shifted to higher values of the disor- 
der. The effect becomes most pronounced for strongly 
correlated disorder, i.e., small a. For an initially insulat- 
ing system with weakly-correlated disorder slightly larger 
than Wc, we find a transition to metallic behavior upon 




FIG. 3: The phase diagram of the 3D Anderson model of 
localization for values of a shown in the legend. For a ^ Qc 
the uncorrelated phase diagram is recovered. For clarity, error 
bars of one standard deviation are shown for every third data 
point only. 



decreasing a. In region IV, the phase boundary moves 
into the metallic phase I in contrast to region III. Thus 
states which are extended for weakly correlated disorder, 
become localized for strong correlation a — > 0. As has 
been argued before jH^, the critical behavior near 
the mobility edges is governed by quantum interference 
as well as tunneling between potential wells. For small 
W and large \E\ i.e., at the band edg e the potential wells 
govern the localization behavior |49| . Obviously, a disor- 
der potential smoothened by correlations will also affect 
these tunneling processes, and away from the band cen- 
ter, this results in more localization. A similar effect has 
been seen for the ID case |2^ . 

The different influence of scale-free disorder on the 
states at the band center and at the band edge is also 
the origin of the difference in the behavior of i^£;(a) and 
vwict) shown in Figs. ^ and |21 The extended Harris cri- 
terion essentially states that the critical behavior of the 
MIT remains unaffected if, as E ^ E^ the divergence 
of the localization length ^ ^ \E — Ec\^'^'^ is stronger 
than the divergence of a typical size \E — i?c|~^^" asso- 
ciated with a correlated disorder potential fluctuation of 
energy \E ~ Ec\. Otherwise, ve = 2/a [Uli^. Let us 
approximate Ec{W) close to Wc as 

Ec{W)^\W -Wcf. (6) 

Then at S = 

^^\E-Ec\-''^ = \Ec\-''^ « \W-Wc\-'"'^, (7) 
so that the critical exponents should be related as 

vw ~ Pi's ■ (8) 



From Fig. we see that /3 -C 1 in the band center 
whereas /? = 1 is possible for \E\ > 0. Furthermore, 
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the above Harris-type argument is based on classical po- 
tential fluctuations, whereas at E = the physics is dom- 
inated by quantum interference with an unchanged uni- 
versal Anderson exponent vyy ?» 1.6 [i^ . 

In summary, we have shown that the extended Harris 
criterion is obeyed varying E at fixed W. The resulting 
exponent ve agrees very well with the predictions of the 
extended Harris criterion. This is the first such demon- 
stration for a fully quantum coherent situation in three 
dimensions to the best of our knowledge. Moreover, we 
show that the extended Harris criterion fails at the band 
center and trace this failure to the different mechanisms 
governing the MIT in the vicinity of the band center and 
outside. We emphasize that such different scaling behav- 
ior at the band center and at the band edges has indeed 
been speculatively discussed for a long time ji^ , al- 
though numerical studies of ue and vw for the uncorre- 
lated case suggest a common value. Here we show that 
suitably long-ranged power-law correlations with a < ac 
give rise to a difference in ve and i^w- 
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